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Abstract 

The 3x3 Maki-Nakagawa-Sakata-Pontecorvo (MNSP) lepton flavor mixing matrix may be slightly 
non-unitary if the three active neutrinos are coupled with sterile neutrinos. We show that it is in 
principle possible to test whether the relation \V^i\'^ + |yg2p + = 1 holds or not in a precision 

reactor antineutrino oscillation experiment, such as the recently proposed Daya Bay II experiment. 
We explore three categories of non- unitary effects on the 3x3 MNSP matrix: 1) the indirect effect in 
the (3+3) flavor mixing scenario where the three heavy sterile neutrinos do not take part in neutrino 
oscillations; 2) the direct effect in the (3+1) scenario where the light sterile neutrino can oscillate 
into the active ones; and 3) the interplay of both of them in the (3+1+2) scenario. We find that 
both the zero-distance effect and flavor mixing factors of different oscillation modes can be used to 
determine or constrain the sum of |Vg]^p, |Vg2p and and its possible deviation from one, and 

the active neutrino mixing angles and 0^^ can be cleanly extracted even in the presence of light 
or heavy sterile neutrinos. Some useful analytical results are obtained for each of the three scenarios. 
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1 Introduction 



The fg, and neutrinos are active in the sense that they take part in the standard weak interactions. 
They are significantly different from their corresponding mass eigenstates z^^! ^2 '^s ^ ^ result of non- 
degenerate neutrino masses and large lepton flavor mixing p]. Whether there exist the sterile neutrinos, 
which do not directly take part in the standard weak interactions, has been an open question in particle 
physics and cosmology. One is motivated to consider such "exotic" particles for several reasons. On the 
theoretical side, the canonical (type-I) seesaw mechanism [2j provides an elegant interpretation of the 
small masses of l'^ (for i = 1, 2, 3) with the help of two or three heavy sterile neutrinos, and the latter can 
even help account for the observed matter-antimatter asymmetry of the Universe via the leptogenesis 
mechanism [3]. On the experimental side, the LSND [5], MiniBooNE and reactor [6] antineutrino 
anomalies can all be explained as the active-sterile antineutrino oscillations in the assumption of one or 
two species of sterile antineutrinos whose masses are below 1 eV |7j. Furthermore, a careful analysis of 
the existing data on the Big Bang nucleosynthesis [8] or the cosmic microwave background anisotropy, 
galaxy clustering and supernovae la [9] seems to favor at least one species of sterile neutrinos at the 
sub-eV mass scale. On the other hand, sufficiently long-lived sterile neutrinos in the keV mass range 
might serve for a good candidate for warm dark matter if they were present in the early Universe [10] • 
That is why the study of sterile neutrinos becomes a popular direction in today's neutrino physics [llj. 

In the presence of small active-sterile neutrino mixing, the conventional 3x3 Maki-Nakagawa-Sakata- 
Pontecorvo (MNSP) lepton flavor mixing matrix [12] is just the submatrix of a (3 + n) x (3 + n) unitary 
matrix V which describes the overall flavor mixing of 3 active neutrinos and n sterile neutrinos in the 
basis where the flavor eigenstates of the charged leptons are identified with their mass eigenstates: 
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Hence the 3x3 MNSP matrix itself must be non-unitary. From the point of view of neutrino oscillations, 
one may classify its possible non- unitary effects into three categories: 

• the indirect non-unitary effect arising from the heavy sterile neutrinos which are kinematically 
forbidden to take part in neutrino oscillations; 

• the direct non-unitary effect caused by the light sterile neutrinos which are able to participate in 
neutrino oscillations; 

• the interplay of the direct and indirect non-unitary effects in a fiavor mixing scenario including 
both light and heavy sterile neutrinos. 

In each of the three cases, no matter how small or how large the mass scale of sterile neutrinos could 
be, the experimental information on the matrix elements of V is essentially different from that in the 
standard case (i.e., the case in which V is exactly a 3 x 3 unitary matrix). Hence testing the unitarity 
of the 3x3 MNSP matrix is experimentally important to constrain the flavor mixing parameters of 
possible new physics and can theoretically shed light on the underlying dynamics responsible for the 
neutrino mass generation and lepton flavor mixing (e.g., the 3x3 MNSP matrix is exactly unitary in the 
type-II [13] seesaw mechanism but non- unitary in the type-I [2] and type-Ill [H] seesaw mechanisms). 
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Following a similar strategy for the precision test of the unitarity of the 3x3 Cabibbo-Kobayashi- 
Maskawa (CKM) quark flavor mixing matrix [1], here we concentrate on a possible experimental test of 
the normalization relation jV^j^P + |V^2p ~^ l^aP = 1 of the 3x3 MNSP matrix in a precision reactor 
experiment which is expected to be able to distinguish between the oscillation modes induced by Am|^ 
and The recently proposed Daya Bay II reactor antineutrino oscillation experiment [15] is just of 

this type, so is the proposed RENO-50 reactor experiment [T6]. At present a very preliminary constraint 
on the sum of \V^i\'^, ^'^d l^gp is [17] 

iKiP + 1^2!^ + = 0-994 ± 0.005 (2) 

at the 90% confidence level, implying that the 3x3 MNSP matrix is allowed to be non-unitary at 
the < 1% level. In this note we shall discuss how to determine |Vg]^p, |Vg2p and |Vg3p via a precision 
measurement of the — )• oscillation and examine whether their sum deviates from one or not. To 
be more specific, we are going to consider three typical scenarios of active-sterile neutrino mixing to 
illustrate possible non-unitary effects on the 3x3 MNSP matrix as listed above: 

• The (3+3) fiavor mixing scenario with three heavy sterile neutrinos which indirectly violate the 
unitarity of the 3x3 MNSP matrix; 

• The (3+1) flavor mixing scenario with a single light sterile neutrino which directly violates the 
unitarity of the 3x3 MNSP matrix; 

• The (3+1+2) flavor mixing scenario in which the light and heavy sterile neutrinos violate the 
unitarity of the 3x3 MNSP matrix directly and indirectly, respectively. 

In each case the sum {V^il"^ + |V"g2p + iKsl^ '^^^ be expressed in terms of the active-sterile neutrino 
mixing angles in a given parametrization of the overall (3 + n) x (3 + n) flavor mixing matrix. Taking 
the parametrization proposed in Ref. [18] for example, we shall show that it is possible to determine the 
active neutrino mixing angles and ^^^3 without any contamination coming from the sterile neutrinos. 
We hope that the points to be addressed in the remaining part of this note will be helpful for the 
next-generation reactor antineutrino oscillation experiments, either to test the standard 3x3 MNSP 
flavor mixing picture or to probe new physics via its possible non-unitary effects. 



2 Non-unitary effects 



In the presence of n species of sterile neutrinos, no matter whether they are very light or very heavy, 
the amplitude of the active — )• z/^ oscillation (for a, f3 = e, /x, r) in vacuum can be expressed as 



^ ^AiW+ ^ />,) • Prop(i.,) • A{i.,W- ^ l^) 



E 



V*i exp 



mfL 



(3) 



in which A{W+ l+v^) = V*J^/{VV^)ac. , Prop(i/i) and A{v^W- ^ ) = V^J {VV^)p,3 describe the 
production of z/„ via the weak charged-current interaction, the propagation of free and the detection 
of Vg via the weak charged-current interaction, respectively |17j — |20j (a schematic diagram is shown in 
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Fig. 1 for illustration). Here m- is the mass of the light (active or sterile) neutrino v^, E denotes the 
neutrino beam energy and L stands for the distance between the source and the detector. It is then 
straightforward to calculate the probability Pif^ — >■ Vp) = We obtain 

= E + 2E [^^« [Va^Vp^Vi^V;,) COS A,^. - Im (V^.y^.F^V;,) sin A,^-] , (4) 

i i<j 

where A- = AmfjL/{2E) with Amfj = mj — mj. One may write out the expression of Vipa ~^ '^p) 
or -P(^'cf ~^ '^p) from Eq. (4) by making the replacement V =^ V* . li V is exactly unitary, then 
{W^)^^ = {VV'^)pij = 1 holds and thus 'P{v^ — i^^) = — )■ Vp) is just the conventional formula of 
— )• z^^ oscillations. Here we only pay interest to the reactor antineutrino oscillations in vacuum, 

ni^e^I^e)= (El^e/y-4E(lK/lK,l^sin2^) • (5) 

Note that the energy-independent term on the right-hand side of Eq. (5) is not equal to one if there are 
heavy sterile antineutrinos which do not take part in the oscillation. This point will be made clear in 
the subsequent discussions, in which a few typical active-sterile neutrino mixing scenarios will be taken 
into account. Given the fact that the typical value of £" is a few MeV and that of L is usually less than 
100 km for a realistic experiment, it is completely unnecessary for us to consider the negligibly small 
terrestrial matter effects on V^Vg — )• I7g). 

2.1 The (3 + 3) flavor mixing scenario 

Let us first consider the (3+3) flavor mixing scenario with three heavy sterile neutrinos which indirectly 
violate the unitarity of the 3x3 MNSP matrix. In this case a full parametrization of the 6x6 flavor 
mixing matrix has been given in Ref. |18) . and the elements in its first row read 
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where c^^ = cosO^j and s^j = sinO-e'^^i with 9- being the flavor mixing angles and being the CP- 
violating phases. In view of the fact that the active-sterile neutrino mixing angles 6i^, O^r^ and are 
at most of 0(0.1) in magnitude [T7], the actual values of the three active neutrino mixing angles 0^2; 

and should be very close to those extracted from current neutrino oscillation data by assuming 
the 3x3 MNSP matrix to be exactly unitary. Eq. (6) leads us to 

\Vei? + \ye2? + iKsl' = C?444 - 1 " {^Ia + S?5 + ^fe) , (7) 

^Note that iyV^)aa a-nd {VV'')^^ in front of P{iy^ — >■ Vp) are the normalization factors and can essentially be canceled 
by the same factors coming from the production of at the source and the detection of at the detector [17], which are 
both governed by the weak charged-current interactions. Hence we focus on V{Va — >■ Up) in this work. 
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implying that the sum of |V^iP, |V'g2p and |V^3p is possible to deviate from one either at the percent 
level or at a much lower level. Is such a small effect really detectable in a precision reactor — )• 
oscillation experiment? 

In the present scenario the three hypothetical heavy sterile neutrinos are kinematically forbidden to 
participate in neutrino oscillations. So Eq. (5) can be simplified to 

^ I7J = lo - 4^ sm2 - 4B sm^ 

+2C sin ^ sin ^ + 8C sin^ — sin^ — , (8 

2E 2E AE AE ' ^ ^ 

where Amgg = Amg^^ — has been used, and 

^0 = (IKil' + IK2l' + |v;3l')' , 
A = (IKil' + IK3l')IK2l', 
B = (IKil' + IK2l')IK3l', 

c = IK2l'IK3l' ■ (9) 

Note that the third oscillation term in Eq. (8) is sensitive to the unknown sign of Am|^ (i-e., Am^i > 
and Am|^ < correspond to the normal and inverted neutrino mass hierarchies, respectively), and it 
might be measurable in a precision reactor antineutrino oscillation experiment in the foreseeable future. 
Some comments on the implications of Eqs. (8) and (9) are in order. 

• The non-unitarity of the 3x3 MNSP matrix, characterized by the small deviation of Iq from 
one, can be tested via the energy-independent zero-distance effect p!9] in such a disappearance 
antineutrino oscillation experiment: 

ViPe ^ ^e) Il=0 =Io = 444 - 1 - 2 (S?4 + S?5 + sfe) ■ (10) 

This effect is in principle measurable with the help of a near detector, although it might in practice 
be indistinguishable from the background which at least includes the uncertainties associated with 
the antineutrino flux. If the energy spectrum of the 17^ — )• 17^ oscillation can be fully established, 
however, it should be possible to determine or constrain the size of Iq. 

• The flavor mixing factors A, B and C are simple functions of the matrix elements |V^iP, |^2p 
and |Vg3p. It is therefore straightforward to obtain 



IKil 



2 _ J{ A-C) (B-C) 

c 



,2 _ j {A-C)C 
l^^2l - \J B-C ' 



Provided A, B and C are all measured to a good degree of accuracy, one may also use Eq. (11) 
to calculate the sum \ V^i\'^ + |V"g2p + iKsP — \/^ examine whether it departs from one or 
not. On the other hand, the value of Iq to be determined in this way can be compared with the 
one to be measured from the zero-distance effect as shown in Eq. (10). 
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Given the concise parametrization in Eq. (6), the active neutrino mixing angles ^"^^ ^13 can 
be determined from the ratios A/Iq, ^ I^q ^ I-^q without any contamination coming from the 
heavy sterile neutrinos. This point is clearly seen as follows: 

^ = j(sin2 20i2cos^^i3 + sin2 0i2sin2 20i3) , 

-'0 ^ 

B 1 , 

- = -sin2 20i3, 

-'0 * 



C 1 

— = -sin2 6'i9sin2 26' 



12 ^^13 



(12) 



Taking 6*12 ^ 34° and 6*13 ~ 9° for example [H], we immediately have A/I^ ~ 0.212, B/Iq ~ 0.024 
and C//g ~ 0.0075. The smallness of C makes it rather difficult to be measured (in other words, a 
determination of the sign of Amg^ or the neutrino mass ordering is a big challenge to the reactor 
antineutrino oscillation experiments [22] — |26j). 



To observe the non-unitary effects in the (3+3) flavor mixing scenario under consideration, one may also 
study the appearance neutrino oscillation (e.g., the — )• oscillation) experiments with reasonably 
long baselines. In such a precision accelerator neutrino oscillation experiment, even new CP-violating 
effects are likely to show up at the 0(10^'^) or 0(10^^) level provided the relevant active-sterile flavor 
mixing angles are not strongly suppressed [27] , 

At this point it is worth mentioning that the type-I seesaw mechanism around the TeV energy 
scale can simply lead to the (3+3) flavor mixing scenario with appreciable non-unitary effects [28j . 
Similar effects may result from some interesting extensions of the type-I seesaw mechanism, such as the 
inverse seesaw scenario [29], the linear seesaw scenario [3^ and the multiple seesaw scenarios [31]. In 
such model-building exercises, presumably significant lepton-flavor-violating effects and non-standard 
neutrino interactions are also expected to show up and their interplay with the non-unitary effects offers 
an important and complementary window for probing the true mechanism of neutrino mass generation 
and lepton flavor mixing |32J. On the experimental side, the precision reactor antineutrino oscillation 
experiments will be complementary to the precision accelerator neutrino experiments for our physics 
goal and thus deserve particular attention. 

2.2 The (3 + 1) flavor mixing scenario 

To illustrate the direct non-unitary effect on the 3x3 MNSP matrix, we consider the (3+1) flavor 
mixing scenario with a single light sterile neutrino motivated by the LSND [4J, MiniBooNE [5] and 
reactor [6] antineutrino anomalies. In this case we have \V^i\'^ + |Vg2p + iKsP — ^~ IK4P — c°s^ by 
using the same parametrization as used above |l8j, where is the active-sterile neutrino mixing angle. 
Because the light sterile antineutrino takes part in the — )• oscillation, Eq. (5) is now simplified to 

ViV^^V^ = 1 - A A' sin^ ^ 4B' sin^ ^ AX' sin^ — 

+2C' sin ^ sin ^ + 8C" sin^ — sin^ — ^ 

2E 2E AE 4E 

+2Y' sin ^ sin ^ + 8Y' sin^ ^ sin^ ^ 

2E 2E 4E AE 

+2Z'sin ^sin ^ + 8Z' sin^ — ^ sin^ — ^ , 13 

2E 2E AE AE ' ^ ' 
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where 

A' = {I- I V;2p) I = ^ (sin^ 2(9i2 cos"^ 6*13 + sin^ 9^^ sin^ '^^13) cos'^ + \ sin^ 6^12 cos^ 6*13 sin^ 26*14 , 

^' = (1 - IK3I') IK3I' = 7 sin2 2^13 cos^ + ] sin^ sin^ 29,, , 
1 



4 — ^•^ ' 4 



= \ye2\>e-f = ^Sm2 0i2Sm2 20i3COS^014 , 

X' = {l-\V,f)\Vj' = ^sm^29,,, 

y = \Ve2\^\Vef = ^sin^ 012 cos^^iasin^ 2^14 , 

Z' = = isin^ sin2 2^14 . (14) 

It is obvious that the oscillation terms associated with 2C' and 2Y' (or 2Z') are sensitive to the signs 
of Am|i and Am|i (or both of them), respectively. This observation makes sense for the place of the 
sterile neutrino in the whole mass spectrum because its mass is in general unnecessary to be larger 
than the mass scale of the three active neutrinos (inparticular if they are nearly degenerate and close 
to 1 eV) . A simple exercise leads us to the solutions □ 

IK2I' = l{i-VT^^ 
1 



IK3I2 = -ii-VT^r^' 



|K4|2 = ^ (1 - yr^Ix^) . (15) 



One may also obtain \V^,\'^ = ^Y' Z' jC and some interesting correlations such as 

c = \(\- \/r^4i^) (1 - ^/r 



45' 



Y' 1 - VI - 4^' 



Z' 1 - VI - 45' 
A few comments on the above results are in order. 



(16) 



Given a light sterile antineutrino whose absolute mass scale is unspecified or Am|i is unknown, 
it is in principle possible to determine or constrain the active-sterile flavor mixing angle ^14 via 
a precision measurement of B' and C' . The deviation of these three parameters from their 
values in the standard case is characterized by nonzero sin^ ^14. Taking ^12 — 34°, ^13 ~ 9° and 
6*14 ~ 10° as a typical example, we obtain ~ 0.208, B' ~ 0.023 and C ~ 0.0071, which can be 
compared with N ~ 0.212, B' ~ 0.024 and C ~ 0.0075 in the 6*14 ~ 0° case. 

Of course, a direct measurement of the oscillation term driven by Am|i will provide the direct 
evidence for the existence of a light sterile antineutrino in the — )■ oscillation. The most 
optimistic case, which seems to be very unlikely, is that X', Y' and Z' could all be determined. 
Considering 6*12 ^ 34°, 6*13 ~ 9° and 9,, ~ 10°, we have X' ~ 0.029, Y' ~ 0.0089 and Z' ~ 0.00072. 
One can see that Z' is too small to be measured, but it does not matter because the sign of Ami-, 



■^Note that the other solution of | V^iP (for i — 2,3, 4) is expected to be inconsistent with the observed neutrino mixing 
pattern, and thus it is ignored here. 



7 



is essentially determinable from a precision measurement of the interference term associated with 
1Y' in Eq. (13). In particular, it should be noted that Y' jZ' = svc? 0^2^'^^ — 12 is a result 
independent of the input value of Q^a,- 

• Provided Aml^ ^ ^"^31 — 2.5 x 10^'^eV^ holds (for i = 1,2,3) and the experimental baseline 
length L favors the oscillation terms driven by Am|]^ and Amgg, the corresponding active-sterile 
antineutrino oscillation terms sin^[Am^^L/(4£')] will practically be averaged to 1/2 because they 
oscillate too quickly. In this case the overall effect induced by the sterile antineutrino becomes the 
zero-distance effect: 

r(pe ^ ^e) Il=o = 1 - 2 (1 - 1^4!') IK4I' = 1 - ^ Sin' 2^14 . (17) 

A near detector is certainly difficult to measure such a small effect because of the uncertainties 
associated with the antineutrino flux. Nevertheless, we remark that the similar effect could be 
cross-checked from a precision determination of the flavor mixing factors A' , B' and C. 

In the same way one may discuss an analogous flavor mixing scenario with two or three light sterile 
antineutrinos which take part in the 17^ — )• oscillation. 

2.3 The (3 + 1 + 2) flavor mixing scenario 

Given three sterile neutrinos, one of them can be assumed to be light enough so as to either interpret 
the existing sub-eV antineutrino anomalies or account for the keV warm dark matter [33]. In such a 
(3+1+2) flavor mixing scenario the two heavy sterile neutrinos may play the role in realizing the seesaw 
and leptogenesis mechanisms [31]. The deviation of the sum |V^xP + |V^2p "I" l^sP from one is also 
described by Eq. (6), but the — )• z7g oscillation probability turns out to be 

V{V^^V^) = /q - 4 A sin^ — 45 sin^ — AX sin^ 



AE AE AE 

+2C sin ^ sin ^ + 8C sin^ ^ sin^ ^ 

AE AE 

+2Y sin ^ sin ^ + sin^ — sin^ 



2E 


sm 


2E 


Aml^L 


sin 


Am\^L 


2E 


2E 


Amg^L 


sin 


Am\^L 



in which 



AE AE 

+2Z sin milMT sin + sin^ ^^u^ -^^2 Am4iL 

2E 2E AE AE ' ^ ^ 



h = (l-lKsl'-lKel')' , 

A= (lKiP + lK3p + lK4l')lK2l' 

B = (|K1|' + IK2P + IK4|')IK3I' 



X = (lKlP + lK2l' + lK3l')l^e 



e4l 



Y = lK2l'lK4l' 



Z = ■ (19) 

We see that Eq. (18) consists of both direct and indirect non- unitary effects on the 3x3 MNSP matrix. 
A straightforward calculation allows us to express | Vgjp (for i = 1, 2, 3, 4) in terms of the flavor mixing 
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factors {A,B,C) and (X,Y,Z) as follows: 



IK3I 





-c 


-Y^ (^B-C-Z^C 


c 




-c 


- r) (i? - (7 - c 


B-C-Z 




-c 


- r) (i? - (7 - c 


A-C-Y 




A- 


C-Y^ (^B-C-Z^C 


1 


[A- 


Y^ (^B-Z^ - C'2 



- ^ ^^ ^ • (20) 

M - rj (b-z \ - C2 

Of course, |V^4p = -JyZ/C holds too|^. Using the neutrino mixing angles given in Eq. (6), we obtain 



/g = cf^cfg as well as 



A 1 1 

— = - (sin^ 2012 cos^ + sin^ 0^2 sin^ 2^]^3) cos^ 0i^ + - s\t? 0^2 cos^ 0i^ sir? 20^^ , 
^0 

B \ 1 

— = - sin^ 20;^3 cos^ 01^ + - sin^ O^^^ s\v? 20^^ , 
Iq 4 

C 1 

— = - sin^ 0^2 sin^ 20^3 cos^ 0i^ , 

£ = ^sin2 2ei4, 
-'0 

Y 1 

— = - sin^ 0^2 cos^ ^13 sin^ '^(^ia > 
-^0 



^ sin^ 6*13 sin^ 2014 > (21) 
Some discussions and remarks are in order. 



Switching off the light sterile neutrino (i.e., setting \V^^\^ = or equivalently X = Y = Z = 
one can easily reproduce Eq. (11) from Eq. (20). If the two heavy sterile neutrinos are switched 
off, then it is straightforward to reproduce Eq. (14) from Eqs. (19) and (21). The interplay 
between the direct non-unitary effect measured by nonzero 0^ and the indirect non-unitary effect 
characterized by nonzero and 0^^ is therefore transparent. In particular, the relationship 

iKiP + IK2P + IK3I' + IK4I' = \l^o = 44 (22) 
holds, and Iq can in principle be determined from the zero-distance effect. 

Provided A/Iq, B/Iq and C//q are all measured in a precision reactor antineutrino experiment, 
it should be possible to determine the flavor mixing angles 0i2, 0i^ and 0^^ in a way independent 



•^Note that the expressions obtained in Eq. (20) may also cover the case of the (3+1) flavor mixing scenario discussed 
in section 2.2. Namely, | V^ij^ can be expressed in terms of {A' , B' , C') and (X', Y' , Z') in the same form as Eq. (20). 
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of the flavor mixing angles ^^^^ and . The same observation is true for the flavor mixing factors 
X /Iq, y /Iq and Z/Iq, although they must be much more difficult to be measured. Of course, the 
situation might be simplified to some extent if the existence of a light sterile antineutrino could be 
established somewhere else (e.g., with the help of a suitable accelerator neutrino (or antineutrino) 
oscillation experiment). But the important point is that the reactor experiment can always provide 
some independent and complementary information on the same physics, no matter whether it is 
new or just conventional. 



In the present flavor mixing scenario, the effective neutrino mass terms of the tritium beta decay 
f H — ;> + e~ + I^g and the neutrinoless double-beta decay A{Z, A) — )- A{Z + 2, N — 2) + 2e~ 
can be expressed as 

1/2 



16 



(23) 



with = 

^i'-'i2''i3 ~^ ^2'^12''13 ~^ '^S'^is being the standard contribution from the three active 

neutrinos and 



("i)ee(Cl4Cl5Cl6) +^4 



'14H5''16 



(24) 



with (m)gg = 'TT'i(C]^2''i3)^ + "^2('^i2''i3)^ + '^3(^i3)^ being the standard contribution from the active 
neutrinos, respectively. Here we have assumed that a possible contribution from the heavy sterile 
neutrinos is negligible [35]. We see that (m)g > (m)g always holds, but it is difficult to judge the 
relative magnitudes of \{m)^^\ and (m)gg because the relevant CP-violating phases may lead to 
more or less (even complete) cancelations of different terms in them [36j . 

If all the three sterile neutrinos are heavy enough, as in the (3+3) flavor mixing scenario, we shall obtain 
the simpler results (m)g = {m) ^cuCi^Ciq and (m)gg = |(m)gg| c\^c\^c\q. 



3 Concluding remarks 

In the era of precision neutrino physics one of the important jobs is to test the unitarity of the 3x3 
MNSP flavor mixing matrix and probe possible new physics which might give rise to some observable 
non- unitary effects on it. Such effects serve as a special example of possible non-standard interactions 
associated with massive neutrinos and their oscillations [57j. Starting from this point of view, we have 
classified three typical categories of non-unitary effects on the MNSP matrix and illustrated them in a 
precision reactor antineutrino oscillation experiment: 1) the indirect effect in the (3+3) flavor mixing 
scenario where the three heavy sterile neutrinos do not participate in neutrino oscillations; 2) the direct 
effect in the (3+1) scenario where the light sterile neutrino can oscillate into the active ones; and 3) the 
interplay of both of them in the (3+1+2) scenario. We have shown that both the zero-distance effect 
and flavor mixing factors of different oscillation modes can be used to determine or constrain the sum 
of |V"g]^p, |Vg2p and |V^e3p possible deviation from one. In addition, the active neutrino mixing 

angles and can be cleanly extracted even in the presence of light or heavy sterile neutrinos. Some 
useful analytical results, which can be applied to a numerical analysis of the future experimental data, 
have been presented for each of the three scenarios under consideration. 
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We expect that some points of view in this work will be helpful for a realistic reactor antineutrino 
oscillation experiment, such as the proposed Daya Bay II experiment which aims to measure the neutrino 
mass ordering and test the unitarity of the 3x3 MNSP matrix. We admit the big challenges that one 
has to face in practice, either in measuring the zero-distance effect or in detecting different oscillation 
modes, or in both of them. Nevertheless, we strongly hope that the Daya Bay II experiment and other 
precision reactor antineutrino oscillation experiments may play an important role in probing possible 
new physics in the Icpton sector, or can at least be complementary to the precision accelerator neutrino 
oscillation experiments in this respect in the foreseeable future. 

An obvious drawback of this work is the lack of a numerical analysis of the future experimental 
sensitivities to the potential non-unitary effects on the MNSP matrix. This is simply because the main 
characteristics of the currently proposed Daya Bay II and RENO-50 experiments remain too preliminary 
and incomplete, and the existing constraints on the non-unitary effects are also preliminary and more 
or less dependent on the hypothetical active-sterile flavor mixing scheme. In this case wc have to focus 
on the analytical discussions about the salient features of direct and indirect non-unitary effects (and 
their interplay) in the present work, leaving a quantitative study of the same topic in a future work. 

The author is indebted to Y.F. Li and S. Luo for intensive and valuable discussions, and in particular 
to S. Luo for her partial involvement at the early stage of this work. He is also grateful to C. Giunti 
and S. Zhou for helpful comments on the interplay of direct and indirect non- unitary effects. This work 
is supported in part by the National Natural Science Foundation of China under Grant No. 11135009. 
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Figure 1: The schematic diagram for the production of the neutrino via the weak charged-current 
interaction, the propagation of free neutrinos in vacuum and the detection of the neutrino via the 
weak charged-current interaction in the — >• z^^ oscihation process (for a, /3 = e,n,T), in which there 
might be possible non-unitary effects on the 3x3 MNSP flavor mixing matrix due to the presence of 
light and (or) heavy sterile neutrinos. 
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